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Abstract 

One of the legacies of the Chandra era is the discovery of AGN-infiated X-ray cavities in virtually all cool- 
core clusters, with mechanical luminosities comparable to or larger than the cluster cooling rate, suggesting 
that AGN might be responsible for heating clusters. This discovery poses a new set of questions that cannot 
be addressed by X-ray imaging or modeling alone: Are AGN actually responsible for halting cooling flows? 
How is the AGN energy transferred to heat? How tight is the observed balance between heating and cooling? 
Using numerical simulations and a new virtual X-ray observatory tool, we demonstrate that high-resolution, 
high-throughput X-ray spectroscopy can address these questions and that the International X-ray Observatory 
IXO will have the necessary capabilities to deliver these measurements. 

Subject headings: techniques:spectroscopic, X-rays:galaxies:clusters, ISM:kinematics and dynamics, galax- 
ies:jes 



1. INTRODUCTION 

One of the surprising legacies of the Chandra era is the 
discovery o f X-ray cavities in virtually every coo l core galaxy 
cluster (e.g. lBirzan et al.ll2004t [Alien et al.ll2006l) . This led to 
the suggestion that AGN heating is responsible for halting the 
cooling in clusters and keeping cooling flows from actually 
reaching star-forming temperatures. 

The discovery of the cavities allowed estimating the kinetic 
power of a large sample of radio galaxies, to within about an 
order of magnitude. The remaining uncertainties result from 
the fact that the cavity age cannot be determined from imag- 
ing studies alone. A range of age estimates have been used 
in the literature, ranging from the sound crossing time to the 
buoyant rise time (e .g. McNamar a et al.l 120 00; Fabia n etldl 
l2000t [Churazo v et al.ll2001l) . Additionally, the orientation of 
the cavities relative to the line-of-sight cannot be determined 
from imaging observations alone, lea ding to an additional un - 
certainty in cavity size and age (e.g. lEnBlin & Heinz! |2002|) . 
Finally, given that we do not know the inflation history of 
cavities, it is not clear how much energy has to be expended 
by the black hole to inflate a cavity to a given size. 

Thus, while it is now generally accepted that AGN can re- 
lease enough energy in principle to heat clusters (at least to 
within a factor of order unity), it is not clear whether and how 
this energy is transferred to the thermal gas in order to actually 
counteract the cooling. 

To make progress, a quantitative spectroscopic approach is 
necessary. Unfortunately, neither Chandra nor XMM-Newton 
have the high spectral resolution and throughput required to 
observe the physical processes that control feedback. As we 
will show below, the requirements to solve the outstanding 
problems in cluster feedback posed by Chandra will be met 
by the International X-ray Observatory, IXO. Before present- 
ing our simulations, it is worth reviewing the relevant techni- 
cal specifications currently planned for IXO: 

• Spectral resolution: IXO will feature an X-ray Micro- 
calorimeter instrument (XMS) with uniform 2.5eV 
spectral resolution up to 7 keV. At the Fe Ka line en- 
ergy of ~ 6.7 keV, this corresponds to a spectral resolv- 
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ing power of 2700, an improvement of roughly two or- 
ders of magnitude compared to the Chandra resolution 
for extended sources. The core XMS array will have 3 
arc-second pixels, arranged in a 40x40 square for a 2x2 
arc-minute field-of-view. 

• Effective area: Relative to Chandra imaging with the 
back-illuminated S3 chip, the effective area will be 
roughly a factor 25 higher, corresponding to the effi- 
ciency difference between the Keck telescopes and the 
Hubble Space Telescope HST. 

• Angular resolution: The angular resolution of the IXO 
telescope mirror assembly will be 5 arc-second on-axis 
at 7 keV (5 arc-second half-power diameter), which is 
very well suited for cluster studies. Again this is com- 
parable to the angular resolution achievable by ground- 
based large telescopes compared to HST 1 . 

The organization of this letter is as follows: In $2] we de- 
scribe the methods employed in the simulations presented in 
$3] while Sj4]discussed the consequences and ^presents our 
conclusions. Throughout this letter, we use concordance cos- 
mological parameters of Hq — 70kms _1 Mpc -3 , tt m = 0.3, 
and tt\ = 0.7 and cluster metallicities of 0.5 relative to solar. 

2. SIMULATIONS 

A powerful way to make detailed observational predictions 
of feedback is by constructing dedicated numerical simula- 
tions of radio galaxies in clusters and then virtually observing 
them with an X-ray telescope. This is the avenue chosen in 
this letter. Before presenting our results, we will briefly de- 
scribe the method used. 

2.1. 77ze code 

We use the publicly available FLASH code dFryxell et al.1 
2000) which is a modular block-structured adaptive mesh re- 
finement code. It solves the Riemann problem on a Carte- 
sian grid using the Piecewise-Parabolic Method. Our simula- 
tion includes 7 x 10 5 dark matter particles. The particles are 

1 While the 0.5 abscond Chandra resolution is necessary for the study of 
very fine details like the sound waves in Perseus, a fair fraction of cluster 
science enabled by Chandra could be done at 5 arc-second resolution. 
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advanced using a cosmological variable-timestep leapfrog- 
method. Gravity is computed by solving Poisson's equation 
with a multigrid method using isolated boundary conditions. 
Our simulations have an effective dynamic range between 
10000 and 20000. 

In order to simulate radio galaxies in clusters, we in- 
ject back-to-back supersonic jets by placing arbitrary inflow- 
boundary conditions into the grid. The location of the "jet 
nozzle" tracks the dynamical center of the cluster potential. 
The jets are injected with a velocity of 0.1c and an internal 
Mach number of 37. 

In order to reproduce the observed morphologies we impose 
a sub-grid jitter on the jet, confined to a 20 degree o pening 
angle, known as the dentist drill effect (Scheuer 1982). More 
information on the numerical detail s of the type of simulations 
used in this letter are presented in Hei nz et all ([2006). 

2.2. Initial conditions 

We used the S2 galaxy cluster from Springel et al. (200j]) as 
our initial conditions. This is a cosmologically evolved rich 
cluster with a mass of M ~ 10 15 M Q . We imported the ini- 
tial conditions into FLASH and evolved the cluster for several 
dynamical times to eliminate any transients before setting off 
the jet. 

The density of the S2 cluster is appropriate for all the clus- 
ters presented in this paper. The temperature profile is ap- 
propriate for our simulations of Cygnus A and Hydra A, but 
we had to adjust the temperature normalization for our simu- 
lations of the Perseus cluster, which has a moderately lower 
temperature. 

Because the cluster is fully cosmologically evolved, our 
simulations incorporate a realistic level of density and tem- 
perature sub-structure, as well as a realistic cluster velocity 
field, including turbulence. We did not, however, include a 
sub-grid turbulence model. 

2.3. XIM: a virtual observatory 

The simulation output was then virtually observed using the 
publi cly available in-house tool XIM (see iHeinz & Br uggen 
2009). Taking input grids from numerical hydrodynamic sim- 
ulations, XIM performs spectral modeling of thermal emis- 
sion, including Doppler shifts and ionization balance. It then 
performs spectral projection along an arbitrary line-of-sight, 
PSF convolution, telescope and detector efficiency, and spec- 
tral convolution with the detector response (using the proper 
response files for current and future telescopes). Finally, it 
adds sky- and instrument backgrounds and Poisson counting 
error. 

2.4. Fe XXV as a kinematic tracer 

In this paper we concentrate on the Ka line from Helium- 
like iron Fe XXV. This line complex is very luminous in the 
intra-cluster gas and because of the large mass of Fe relative 
to other abundant sp ecies, it is an excellent kinematic tracer 
dBriiggen et alj|2005l) . 

Many other lines are present, also showing the kinematic 
signatures discussed below, but for simplicity and clarity we 
limited the analysis to Fe XXV. Figure[T]shows the rest-frame 
multiplet line structure at a plasma temperature of 4 keV for 
reference. Given its strength, the resonance line at 6.7 keV 
rest energy is the best suited for kinematic measurements. 

3. RESULTS 
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FIG. 1 . — Reference spectrum of the Fe XXV Ka line for a 4keV thermal 
plasma, with resonance, intercombination, and forbidden lines labeled. Di- 
electronic satellite lines are also included. 

To investigate the possibilities offered by high-throughput 
X-ray spectrographs for the study of AGN feedback in clus- 
ters, we constructed virtual IXO observations of three repre- 
sentative (benchmark) galaxy clusters with detailed Chandra 
cavity studies: Perseus A, Cygnus A, and Hydra A. The sim- 
ulations aim to reconstruct, as best as currently possible, the 
observational characteristic of these sources. 

Before presenting the three cases in more detail, it is 
worth pointing out that IXO data are difficult to visualize 
in 2D figures. The plots show different spectra from vir- 
tual spectral slits placed across that data cube, which clearly 
show the ability of IXO to resolve the kinematic signa- 
tures of cluster feedback. However, a much better way to 
view the data is by animation. For this reason we have 
prepared movies of all three cases that can be viewed at 
http://www.astro. wise. edu/^heinzs/feedbackmovies. html 

3.1. Perseus A 

The Perseus cluster was the first cluster known to have X- 
ray cavities. With a mega-second of Chandra observing time, 
this is the best studied cluster to date, displayin g rich morpho- 
logical complexity across all wavebands (e .g. iGraham et alJ 
l2008tlFabian et al.ll2008t iTavlor et al.ll2006l) . As such, it rep- 
resents the benchmark for numerical models. 

As a moderately powerful source, with a jet power some- 
where between 10 44 and 10 45 ergss _1 , it also represents a 
good example of what ha s become known as "ge ntle" or "ef- 
fervescent" heating (e.g. Ruszkowski et al. 2004). 

The detection of a weak shock in the long Chandra image 
dGraham et al.l l2008) has provided the most robust power esti- 
mate of Wjct ~ 10 45 ergs s , which is the number we chose 
to use. Given this power estimate, we ran the simulation for 
10 7 yrs until the size of the shock and cavities corresponded 
to the observed values. 

The central cluster temperature is around 3-4 keV, some- 
what colder than the S2 cluster we used as our initial con- 
dition. We thus adjusted the cluster temperature in post- 
processing to ensure the correct line strengths and ratios, fix- 
ing the total X-ray power from the cluster to the observed 
Chandra flux. 

The resulting simulated IXO image and Fe XXV Ka spec- 
trum are shown in Fig.|2]for an assumed exposure time of 250 
ksec, at a red-shift of z = 0.01756. The cavities are clearly 
visible, as is the classical bright central core of the Perseus 
cluster. 

Figure [2] also show the spectrum for a virtual spectral slit 
placed across both cavities (see image for slit placement). 
Apart from the complex sub-structure of the Fe XXV Ka mul- 
tiplet, the most striking feature immediately visible from the 
spectrum is how the lines splits at the locations of the cavities. 

The two "bubbles" seen in the spectrum correspond to the 
front and back walls of the cavities, allowing full 3D recon- 
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FIG. 2. — Virtual 250 ksec IXO observation of Perseus cluster (redshift 
z = 0.01756). Top: 0.5-10 keV image. Bottom: Fe XXV Ka spectrum 
of a 2-pixel vertical slit across the cavities, clearly showing the kinematic 
signature of the expanding bubble . 



struction from the imaging spectrum like in the case of su- 
pernova remnants. Simply reading off the expansion velocity 
from the blue- and red-shifted line gives a line-of-sight veloc- 
ity of i>los ~ ±350kms _1 , compared to the actual velocity 
of ~ 325 kms -1 in the simulation. 

Thus, IXO will be able to clearly resolve the kinematic sig- 
nature of expanding bubbles around cluster radio sources like 
Perseus A. It is also clear that angular resolution of a few arc- 
seconds is vital to resolve the cavities both spectrally and spa- 
tially. 

3.2. Cygnus A 

As a powerful FR2 source, Cygnus A sits at the upper 
end of the expected range of cluster radio sources and shows 
clear cavities and p ossibly a shock in the Chandra image 
dWilson et al.l l2006). As with most cluster sources, the am- 
biguities in the interpretation of the imaging data cannot con- 
strain the jet power to better than about an order of magnitude, 
roughly 10 45 to 10 46 ergs s _1 . 

The discovery of giant X-ray cavities in other clusters has 
shown that powerful outbursts like the one currently observed 
in Cygnus A might be more common than implied by models 
of "gentle" cluster heating, making this source an important 
benchmark. 

Our initia l numerical simula tions of Cygnus A have been 
presented in lHeinz et al.l (120061) . Here we will present the vir- 
tual IXO observation derived from that simulation. We ran the 
simulation with a jet power of 10 46 ergs s _1 for 21 Myrs, at 



40 ; 



30 : 



20 : 



10 ; 



oL 





* 



10 



20 



30 40 
ra [pixel] 



50 



60 




6.26 6.28 



6.30 



6.32 
E [keV] 



6.34 



6.36 



6.38 



which point the simulation reached the observed cavity and 



FIG. 3. — Virtual 250 ksec IXO observation of Cygnus A (redshift z = 
0.056). Top: 0.5-10 keV image. Middle: Fe XXV Ka line of the eastern 
cavity. Bottom: same for the western cavity. Both cavities will be easily 
resolved by IXO . 

radio lobe size and morphology. 2 

Our virtual 250 ksec IXO observation is shown in Fig. [3] for 
a cluster redshift of z = 0.056. The two bottom spectra show 
the FeXXV K-alpha line, clearly resolving the two cavities 
(the two spectra correspond to two virtual spectral slits, one 
for each cavity). As expected, for a source this powerful it 
will be easy for IXO to resolve the kinematic structure of the 
line. Thus, IXO will be able to conclusively determine the jet 
power of Cygnus A. 

3.3. Hydra A 

Several clusters that currently contain relatively modest ra- 
dio sources show fossil cavities that imply much more pow - 
erful past o utbursts, most no tably Hydra A dWise et aLll2007f). 
MS 735 dGitti et al.l l2007h . and Hercules A dNulsen et alJ 
120051) . Confirming this interpretation will be an important 
test of impulsive vs. gentle cluster heating models (note that 
IWise et aD (120071) identify multiple powerful outbursts in Hy- 
dra A). 

We used a late timestep of our 10 46 ergss _1 simulation 
to represent this case, 150 Myrs after the outburst started, 
roughly corresponding to the case of Hydra A. The cavity 
sizes are of the same order as those observed in Hydra A, in- 
cluding the north-south asymmetry due to motion of the intra- 
cluster medium (ICM). 

Given the low surface brightness this far out in the cluster 
and the large angular size of these sources, such observations 
will be more challenging (requiring multiple long pointings). 
Still, the 250 ksec virtual observation in Fig.@]shows that IXO 
can resolve and detect the velocity structure of the cavities in 
these sources. 

2 We also tested the predicted virtual Chandra surface brightness of the 
simulation against the actual observation and found that they agree to within 
better than 10% in the region of the rim around the cavity, giving us confi- 
dence in the predicted line fluxes. 
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FIG. 4. — Virtual 250 ksec IXO observation of a Hydra A-like radio source. 
Top: 0.5-10 keV image. Middle and bottom: Fe XXV Ka spectra of both 
cavities. 4. DISCUSSION 

4.1. Implications for other cluster radio sources 

The three detailed examples presented above show that a 
IXO will be able to resolve the kinematic structure of radio 
galaxy drive cavities for a significant fraction of the sources. 
In the case of Perseus A, the large count rate will allow a 
clear resolution of the iron line in observations as short as 20 
ksec. This implies that even in clusters with lower surface 
brightness, IXO would be able to deliver cavity velocities in 
exposures of moderate length (100-200 ksec). 

In principle, IXO will be able to resolve cavities expanding 
with velocities as slow as i>los <L 100 km s -1 , though the 
presence of turbulent broadening in the ICM will likely mean 
a higher threshold (see ij4.2t . Given that the expected expan- 
sion velocity of actively jet- driven cavities for a fixed physical 
cavity radius R scales like dHeinz et alj [l998) 

/w^ 1/3 
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low-power sources like M87 would be more difficult to ob- 
serve. Detailed predictions of what IXO would observe in the 
case of M87 will have to await more dedicated simulations. 

4.2. Caveats 

Given that the presented virtual observations are based on 
numerical simulations, there is are several physical effects 
that we did not include and that could in principal affect the 
predicted results. However, the point of this paper is not the 
precise prediction of velocities of individual sources but rather 
to demonstrate that IXO will be able to measure the expansion 
rates for the expected range of velocities. We will therefore 
not go into any detail of the numerical shortcomings (like the 
absence of magnetic fields, the limitation to non-relativistic 
simulations, the restriction to purely thermal particle popula- 
tions, and neglecting effects of radiative transfer like resonant 
scattering). 



One of the critical elements that should be addressed is the 
level of turbulence present in the ICM significant turbulence 
will broaden all lines and make kinem atic measurements of 
the sort discussed here mor e difficult (flnogamov & Sunvaevl 
12003b IChurazov etai1l2004l) . While our simulations do in- 
clude turbulence on the scales resolved by our grid, we can- 
not make any statements regarding the level of sub-grid tur- 
bulence. Turbulent velocities in excess of 300 km s -1 will 
impede cavity measurements. 

Finally, calibration uncertainties present in any X-ray in- 
strument will in principle affect the spectroscopic accuracy. 
However, because we are only discussing kinematic effects, 
systematic errors due to calibration errors should not have any 
strong effect on our results. 

4.3. Other prospects for feedback studies 

The analysis presented above is only the most forward and 
simple application of high-resolution, high-throughput X-ray 
spectroscopy to the questions of cluster feedback. Given the 
spectral complexity of typical thermal spectra, many more ap- 
plications will be possible that are well beyond the scope of 
this letter. We will suggest a few avenues that promise to be 
particularly fruitful: 

• Line stacking: Given IXO 's high spectral resolution 
over a wide bandwidth, it will be easy to re-grid the 
spectral data to a logarithmic grid, upon which one can 
stack different emission lines on top of each other sim- 
ply by shifting along the energy axis, greatly increas- 
ing signal-to-noise (thus reducing the overall exposure 
time requirements). More complex deconvolution al- 
gorithms could yield even more powerful diagnostics, 
including the temperature structure of the cavity shells 
from line-ratio variations. 

• Line-of-sight angles: From the velocity centroid of the 
two cavities, relative to the cluster mean, it will be pos- 
sible to determine the mean jet orientation relative to 
the line-of-sight as well as detecting large scale rota- 
tion in the cluster. 

• Emission measure mapping: Beyond simple kine- 
matics, spectral fitting will allow the extraction of the 
emission measure distribution in each individual X-ray 
pixel, allowing multi-phase metallicity, entropy, tem- 
perature, and pressure maps. The ultimate promise of 
this approach is the direct detection of the spectroscopic 
signature of AGN heating. The performance require- 
ments on the XMS calibration to allow such a measure- 
ment has yet to be formulated. 

5. CONCLUSIONS 

We showed that the quantitative kinematic measurements 
critical for establishing accurate cavity ages, jet powers, and 
black hole duty cycles are only achievable through high- 
resolution X-ray spectroscopy. Our simulations demonstrate 
that the 2.5 eV spectral resolution and 5 arc-second spatial 
resolution provided by the International X-ray Observatory 
IXO will be sufficient to make the necessary measurements. 
Any significant compromise in angular or spectral resolu- 
tion would render these measurements impossible, implying 
that neither the current or any other planned X-ray telescopes 
(such as Astro-H) could deliver them instead of IXO . 
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